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Abstract 



. The complete 0{a) corrections including soft-photon bremsstrahlung to the pro- 

J> I cess e~^e~ — > W~^W~ in the MSSM are calculated for on-shell W bosons. The 

relative difference between the MSSM and Standard Model corrections is generally 
. quite small. The maximum deviation from the Standard Model within the scanned 

i region of parameter space is < 1.5% for unpolarized and transversally polarized W 

^ I bosons, and < 2.7% for longitudinal W bosons. 

O ; PACS numbers: 12.60.Jv, 13.10. +q, 12.15.Lk. 

Oh! 1 Introduction 

The process e~^e W^W is already one of the key processes at LEP2, and will be 
of similar importance at future linear e^e~ colliders. Hence it is not surprising that con- 
siderable theoretical effort has gone into the precise prediction of the cross-section in the 
5^ ! Standard Model (SM), both for on- and off-shell W bosons see for a review). 

For a process well accessible both experimentally and theoretically in the SM, one of 
the obvious questions to ask is whether it can tell us anything about physics beyond the 
SM. Super symmetric extensions play a special role because they, like the SM, allow to 
make precise predictions in terms of a set of input parameters. Previous calculations in 
supersymmetric theories include the complete one-loop corrections in spontaneously broken 
supersymmetry 0, sfermion-loop effects in the MSSM 0, ^, and also the complete MSSM 
corrections to the closely related triple-gauge-boson vertex 

In this paper the complete one-loop corrections for e^e^ W^W~ in the MSSM 
including real bremsstrahlung in the soft-photon approximation are presented. The results 
are rather small: the maximum deviation from the Standard Model within the scanned 
region of parameter space is < 1.5% for unpolarized and transversally polarized W bosons, 
and < 2.7% for longitudinal W bosons. 

The outline of this paper is as follows. In Sect. ^ the kinematics and notation are fixed. 
Sect. ^ describes the details of the one- loop calculation. The results of the calculation and 
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the scan over the MSSM parameter space are presented in Sect. ^. Sect, ^finally gives the 
conclusions. 



2 Kinematics and notation 

The reaction studied here is 

e+(fci, Ai) + e-(fc2, A2) ^ W^{h, A3) + W-{k^, A4) , (1) 

where fcj and Aj represent the momenta and hehcities of the external particles, respectively. 

The incoming particles travel along the z axis and are scattered into the x-z plane. Ne- 
glecting the electron mass, the explicit representations of the momenta and the polarization 
vectors in the centre-of-mass system are 

= E{1, 0, 0, -P) , 4(0) = {-p, EsinO, 0, E cos 9) /Mw, 

k^ = E{1, 0, 0, (3), 4(±) = (0, -cos^, ±i, sm9)/V2, 

k^ = {E, -p sin 9, 0, -p cos 9) , e^(0) = {-p, ~E sin 9, 0, -E cos 9) /Mw , 



^4 = {E, p sin 9, 0, p cos 9) , ^ti^) — (0' '^'^s ~ ^in 6 

where P is the velocity of the electrons, p = E"^ — is the momentum of the W bosons, 
E = 1/5/2 is the beam energy, and 9 is the scattering angle. 

The polarized differential cross-section is obtained from the helicity amplitudes A4 as 

/ da A -'- I A ^ |2 

In the following, only unpolarized electrons are considered, since that is the situation at 
LEP2. In this case the cross-section has to be averaged over the initial helicities according 
to 

^) ^ (I) ■ w 

'^"/UUA3A4 ^ XiM=±l ^AiA2A3A4 

Polarization combinations are given by a sequence of four letters, e.g. UULT, where L, T, 
and U denote longitudinal (A = 0), transverse (A = ±), and unpolarized particles. 



3 One-loop corrections 

Since the coupling of the electron to any of the Higgs particles in the SM or the MSSM is 
suppressed by a factor rrie/Mw, the Higgs-exchange diagrams can safely be neglected. The 
tree-level diagrams are then the same for the SM and the MSSM, i.e. 7 and Z exchange 
in the s-channel and neutrino exchange in the t-channel (see Fig. |^). The tree- level results 
shall not be discussed here as this has been done in detail elsewhere (explicit formulas can 
be found e.g. in |0]). 
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Figure 1: The tree- level diagrams. 
3.1 Calculational framework 

The calculation of the 0{a) radiative corrections has been performed in 't Hooft-Feynman 
gauge. Ultraviolet (UV) divergences were treated within dimensional regularization. For 
the renormalization the on-shell scheme was used, following the formulation worked out 
in p. In the absence of SUSY particles at tree level, only those counter-terms appear that 
are already present in the SM. The only change is that now the self-energies from which 
the renormalization constants are derived have to be calculated in the MSSM. 
To 0{a), the squared matrix element is given by 

\M\^ = |^Bornr(l + ^soft) + 2 Rc (A<^„,,A<i.loop) , (5) 

where A^Bom and A^i-ioop denote the sum of the contributing tree-level and one- loop Feyn- 
man diagrams, respectively, and (5soft is the QED correction factor from real bremsstrahlung 
in the soft-photon approximation. 

The Feynman diagrams were generated with FeynArts which in its current version 
uses algorithms that can deal with supersymmetric theories . The resulting amplitudes 



were algebraically simplified using FormCalc [11 1 and then converted to a Fortran pro- 



gram. The LoopTools package [|TT|, 12 1 was used to evaluate the one- loop scalar and tensor 
integrals. For a single point in phase and parameter space the Fortran program runs about 
3.6 ms in the SM and 102 ms in the MSSM. Scans over parameter space in the MSSM 
typically take several hours. 

Apart from obvious checks such as UV- and IR-finiteness, the SM results of were 
fully confirmed and the effects of sfermion-loops in g| could be reproduced for various 
scenarios of sfermion masses to within a very small constant shift (~ .5%) which is an 
effect of the different renormalization schemes used and hence of higher order. 

With the presence of external gauge bosons one might expect sizable gauge cancella- 
tions and hence instabilities in the numerical code. To safeguard against such problems, the 
SM values were computed once in the conventional way and once in the background-field 
method [|l3l, where gauge cancellations should be greatly reduced, and an agreement to 9 
digits was found. Even though the same comparison cannot easily be done for the MSSM 
(a corresponding model file is currently not available), a similar stability is expected since 
the gauge structure of the MSSM is essentially the same as in the SM. As an alternative 
measure of stability, one can compare the individual contributions from self-energies, ver- 
tices, and boxes with their sum. To show that the cancellations are harmless numerically. 
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the following table lists the individual contributions (each of them renormalized) that make 
up the amplitude at a particular point in parameter space for longitudinal W bosons, for 
which the largest cancellations can be expected. 



647r2, 



v^=500 GcV 
0=90° 



X = self-energies: 
X = vertices: 
X = boxes: 



-.02535074 pb 
-.03675248 pb 
.06422023 pb 



X 



sum: 



.00211700 pb 



(6) 



3.2 QED corrections 

In addition to the virtual diagrams, real photon emission from the external legs has to be 
taken into account to cancel the infrared (IR) divergences which arise from the exchange 
of massless photons. The IR divergences are regularized by an infinitesimal photon mass 
A. In the soft-photon limit the cross-section for real photon emission is proportional to the 
Born cross-section, 

=(^) ^so., (7) 



with the soft-photon factor Sso{t given by 



Jsoft 



(27r)3 J 2ko ^ {kik){kjk) 



Here AE is the maximum energy of the emitted photons, Qi is the charge of the ith external 
particle, and the sign in front of the product QiQj is -|- if particles i and j are both either 
incoming or outgoing, and — otherwise. The basic integrals needed for the evaluation of 
(|) have been worked out e.g. in [|I^. 

The absolute magnitude of the 0{a) corrections is largely determined by the QED 
contributions and depends strongly on the photon-energy cutoff AE through logarithms 
of the form \ogAE/E. Without an experimentally motivated choice of AE, the overall 
size of the one-loop corrections can be shifted more or less at will using different values 
of AE. Since the intention here is only to show the differences between the SM and the 
MSSM calculations, the concrete value of AE is rather unimportant. In this calculation a 
soft-photon cutoff energy of AE = 0.05^^5 has been used as in 0. 

The cancellation of the IR divergences has been checked numerically by varying the 
photon-mass A on which the result must not depend. A variation of A from 1 to 10^*^ GeV 
leaves the first 10 digits of the cross-section invariant. 

While soft-photon radiation is sufficient to cancel IR divergences, it is an approximation, 
valid only if AE is small compared to all relevant energy scales. If this is not the case, hard- 
photon radiation must be taken into account, too. However, the hard-photon corrections 
at 0{a) are exactly the same as in the SM, therefore they have been omitted in the present 
calculation. 
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Figure 2: The MSSM self-energy corrections. Braces indicate that there is one diagram 
for the first members of all braced lists, one for the second members, etc. A sfermion with 
index i accounts for six particles, e.g. Cj = {e^, e^, /i^, /i^, f \ f ^}. 

3.3 Inventory of one-loop diagrams 

This section lists the contributing one-loop diagrams. The diagrams can be separated 
into self-energy contributions, vertex corrections, and box corrections. Note that diagrams 
which are already present in the SM are omitted in the figures. 

The self-energy contributions fall into two categories, corrections to the 7 and Z prop- 
agator in the s-channel, and corrections to the u propagator in the t-channel. These 
diagrams are shown in Fig. |[ 

The vertex diagrams can be grouped into corrections to the initial- and final-state 
vertex in the s-channel (Figs. |^ and and corrections to the t-channel vertices (Fig. |^). 

Finally, the box diagrams are displayed in Fig. ^ 



4 Numerical results 
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Figure 3: The MSSM contributions to the initial-state vertex in the s-channel. 



4.1 Input parameters 

4.1.1 Standard Model parameters 

For the SM parameters the following numerical values are used: 

Mz = 91.1867 GeV, Mw = 80.39 GeV, 



a 

m 
m 



rrir 



137.0359895, 

0.51099907 MeV, m„ = 53.8 MeV, = 53.8 MeV, 

105.658389 MeV, = 1.50 GeV, m, = 150 MeV, 

1777 MeV, rrit = 174 GeV, rrib = 4.7 GeV. 



(9) 



The masses of the up and down quarks are effective parameters which are adjusted such 
that the five-flavour hadronic contribution to Aa is 0.02778 HTSll, i.e. 

„(5) 



Ml) 



a 



J2 



f=u,c,d,s,b 



m 



f 



0.02778 . 



4.1.2 MSSM parameters 

Higgs sector The neutral Higgs sector is fixed by choosing a value for tan/5 and for 
the mass M^o of the CP-odd neutral Higgs boson A^. For the other neutral-Higgs masses, 
which receive significant radiative corrections [0, the two-loop approximation formula of 
1% is used. 

There are only small radiative corrections for the charged Higgs masses and the following 
equation holds for M^o ~ 0{Mw) fPi 



M2 



5aM2 MsusY 



In- 



w 



+ 



3a 



47rs2 A/f2 



sin^P cos^/5 



w 



m 



t.r 



+ 



mr 



sin^/? cos'^P 



2 

3" 



j^MsUSY^ (10) 

'mt,r 



where Msusy is a universal soft-SUSY-breaking mass introduced in the next paragraph 
and rrit^r = mtll + 4as(mj)/37r)~-'^ is the running top mass. 
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Figure 4: The MSSM contributions to the s-channel final-state vertex. Braces indicate that 
there is one diagram for the first members of all braced lists, one for the second members, 
etc. A sfermion with index i accounts for six particles, e.g. Cj = {e^, e^, jl^, fl'^, f^,f^}. 
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Sfermions For simplicity, all soft-SUSY-breaking parameters are assumed equal and 
mixing between sfermion generations is neglected, so that 

M% = Ml = M? = M? = Ml 
Then, the sfermion mass matrix is given by |]T^ ^] 



Au = A^l , 



L 

Al 



E 

Ad 



'111 



MlusY + M| cos(2/3)(// - Qfsl) + 
"^/(^{n,d} - /i{cot^,tan/3}) 



m 



where the elements in braces apply to L 



f 



f m/(v4{„_rf} - /i{cot/5,tan/3}) 
^susY + Ml cos{2p)Qfsl, + m) 

+1 and — |, respectively. 



(12) 



Charginos and Neutralinos The chargino mass matrix |]T9|, ^| 

^_ f Ma y2Mii/sin/? 

and the neutralino mass matrix UTOl ^ 
/ 

Y = 



(13) 



Ml 
M2 

—MzSw cos P MzCw cos P 

y MzSw sin P —MzCw sin P —fx 

are diagonalized with unitary matrices U, V, and such that 

U*XV-^ = diag(m^^,m^J , 
N*YN~^ = diag(mj^o, . . . , m^o) . 



-Mz-sw cos P Mz-sw sin P\ 
MzCw cos P —MzCw sin P 







(14) 



/ 



(15) 



"xp • • • ^ ""xi' 

This diagonalization is done numerically using the subroutines of the LAPACK library 



21 1 . The U{1) gaugino-mass parameter Mi which appears as a further input parameter in 

(16) 



Y is fixed, as usual, by the SUSY-GUT relation 

Ml 



-4^M2 



4.1.3 Parameter scan 

The remaining input parameters are then scanned over the following regions: 



SM 


MSSM 




tan/3 = 


1.5,5,50 




Mao = 


100 .. . 1000 GeV 


Mh = 100 . . . 300 GeV 


fi = 
M2 = 


-1000... 1000 GeV 
100 .. . 1000 GeV 




^SUSY = 


100 .. . 1000 GeV 




Au = Ad = 


^SUSY, 2MsuSY 



(17) 
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When the particle masses are calculated from these input parameters, the Fortran program 
checks whether they are consistent with the current exclusion limits and automatically 
omits already-excluded points in parameter space from the scan. The following bounds are 
used: 



Mf > 80 GeV 22 



Ml ^ 70 GeV 



M-^^l-^ ^ 150 GeV [|2 



M, > 70 GeV 23 



Myfi ^ 85 GeV |24| 



m^o ^ 30 GeV 



^ 90 GeV P 



(1^ 



4.2 Results 

In Fig. the absolute values of the total and differential cross-section in the SM and 
MSSM is shown for three different polarizations: UUUU (all particles unpolarized), UUTT 
(transverse W bosons), and UULL (longitudinal W bosons). The thickness of the curves 
reflects the range of the allowed values within the parameter scan (|I^. 

Fig. P shows the relative corrections in the SM and MSSM with respect to the Born 
cross-section. The overall magnitude of the corrections is largely determined by the well- 



understood QED logarithms, as discussed in Sect. 3.2. The bands show the variation 



of the cross-section within the parameter scan ([T7|) . The bands for the SM and MSSM 
overlap slightly, the MSSM band being the lower one. The two bands are of comparable 
width of the order of a few percent for all polarizations. The dominating polarizations over 
the entire energy range are UU±=F, and this explains why the plots for the UUUU and 
UUTT polarizations in Fig. ^ are quite similar. On the other hand it is not possible to 
trace the origins of the radiative corrections by disentangling the self-energy, vertex, and 
box contributions since in the presence of external gauge bosons there are delicate gauge 
cancellations. 

The relative deviation between the SM central value and the MSSM band is plotted in 
Fig. ^ While the largest corrections (~ —2.7% at 1 TeV) are seen for purely longitudinal W 
bosons, one has to keep in mind that the cross-section for longitudinally polarized W bosons 
is much smaller than for the transverse polarizations. In the transverse polarizations and 
in the unpolarized case the maximum deviation between the SM and the MSSM is roughly 
1.5%. The maximum deviation is reached in all polarizations for light SUSY particles, for 
example in the unpolarized case (UUUU) the lightest higgs, stop, chargino, and neutralino 
masses at the point of maximum deviation are Myfi = 91 GeV, = 165 GeV, = 
117 GeV, and m^o = 52 GeV, respectively. 

The variation of the cross-section with the scanned MSSM parameters is shown in 



Fig. |I0[ The largest variation is connected with the soft-SUSY-breaking mass MgusY and 
confirms the idea of the previous calculations Q that important contributions come from 
the sfermion sector. Nevertheless, the contributions from the other sectors are similar in 
size and cannot be neglected. 

The calculation was checked for a possible scheme dependence by repeating it using the 
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a/ph 




200 400 600 800 1000 200 400 600 800 1000 



Vs/GeV Vs/GeV 




tt/2 t 7r/2 T 

e 6 



Figure 7: The total cross-section and the differential cross-section at ^/s = 200 GeV for the 
SM and the MSSM. The thickness of the curves give the range of allowed values within the 
parameter scan (|1^). The polarizations are denoted by: UUUU (all particles unpolarized) , 
UUTT (transverse W bosons), and UULL (longitudinal W bosons). 
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"■l-loop — C"Bom 
C"Bom 



-10% - 



-15% 




uuuu 



1000 



Vs/GeV 



-20% - 



-25% L 



-10% - 



-15% 




UUTT 



1000 



Vs/GeV 



-20% - 



-25% L 



-10% 



-15% 



-20% - 



-25^ 




1000 



Vs/GeV 



Figure 8: The relative difference between the one-loop corrected and the Born cross-section. 
The cross-hatched bands indicate the minimum and maximum reached within the scan 
over the parameter space (|T7|) . The polarizations are denoted by: UUUU (all particles 
unpolarized) , UUTT (transverse W bosons), and UULL (longitudinal W bosons). 
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Fi gure 9: The relative difference between the SM central value Ci.ioop 

and the MSSM 

bands of Fig. ||. The polarizations are denoted by: UUUU (all particles unpolarized) , 
UUTT (transverse W bosons), and UULL (longitudinal W bosons). 
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Figure 10: The variation of the MSSM cross-section for unpolarized particles with the 
scanned MSSM parameters. The variation of a parameter X is computed as the difference 
between maximum and minimum of the cross-section for fixed X, averaged over all values 
of X. 



input values {a^Gp, Mz} instead of {a, Mz}. This is not entirely straightforward 
since Mw enters not only in loop corrections, but more directly in the kinematics of the 
process. For instance, the threshold of 2Mw is different in the SM and MSSM so that it is 
not admissible to plot the ratio as in Fig. 0. To compute Mvy(a, Gp, Mz), the fit formula 



of [26|, approximating the two-loop result, together with the MSSM contributions to Ap 



| P7| have been used. The results are shown for longitudinal W bosons, for which the MSSM 
corrections are largest, in Fig. |Tl| and differ only slightly from the former result. 



5 Conclusions 

The MSSM corrections to e~^e~ W^W~ are probably too small to be detected at LEP2, 
taking into account that the convolution with decay amplitudes in the final and photon 
radiation in the initial state tends to smear the corrections. 

With the projected accuracy of a linear collider, however, corrections of this magnitude 
will play a role. Yet even a linear collider may have a hard time disentangling the MSSM 
contributions because the one-loop electroweak corrections include Sudakov logarithms of 
the form log (s/m^), whereas for SUSY corrections these logarithms are power-suppressed 
28| . Extrapolating from 20% one-loop electroweak corrections at TeV energies, the higher- 



order electroweak contributions may be comparable in magnitude to the one-loop MSSM 
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C"l-loop — CTB orn 
CBorn 



200 



400 



600 



300 



1000 



Vs/GeV 



-10% 



-15% 



-20% - 



-25% 




Figure 11: The same as the lower diagram in Fig. |^ (i.e. the relative one-loop corrections for 
longitudinal W bosons), but calculated from the input parameter set {a, Gp, Mz} instead 
of {a, Mw, Mz}- The cross-hatched bands indicate the minimum and maximum reached 
within the scan over the parameter space (0). 



corrections. 

Although the calculation presented here is limited to on-shell W bosons and hard- 
bremsstrahlung effects have not been considered, these effects are the same for the SM and 
the MSSM, hence it should be possible to use the MSSM cross-section presented here in 
a straightforward way in existing Monte Carlo generators, for example [^. The Fortran 
code for this process is available at [tittp : //www . hep-processes . de| . 
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